SUBROUTINE GRUNG1 DOUBLE PRECISION TD,ZO,EMA,PI,W,WC,X,Y DIMENSION X(50),Y(50) COMMON /GRUNGE/I,TD,ZO,EMA,PI,W,WC COMMON /FIT/ X,Y X(I)=ZO*ZO*W*W Y(I)=ZO*ZO RETURN END